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Low temperature scanning tunneling microscopy is used to image the doped topological insulator 
Fe x Bi 2 -xSe 3 . Interstitial Fe defects allow the detection of quasiparticle interference (QPI), and 
the reconstruction of the empty state band structure. Quantitative comparison between measured 
data and density functional theory calculations reveals the unexpected coexistence of quantum well 
states (QWS) with topological surface states (TSS) on the atomically clean surface of Fe cc Bi 2 - cc Se 3 . 
Spectroscopic measurements quantify the breakdown of linear dispersion due to hexagonal warping. 
Nonetheless, both QWS and TSS remain spin-polarized and protected from backscattering to almost 
1 eV above the Dirac point, suggesting their utility for spin-based applications. 

PACS numbers: 68.37.Ef, 71.20.-b, 73.20.At, 71.15.Mb 


Topological insulators (TIs), recently discovered ma¬ 
terials with insulating bulk and topologically protected 
helical Dirac surface states, have generated widespread 
excitement due to proposed applications such as dissi¬ 
pationless spintronics, ambipolar transistors, and fault- 
tolerant quantum computers [1-3]. Many of these appli¬ 
cations hinge on the predicted ability to preserve spin in¬ 
formation without backscattering. Bi 2 Se 3 has attracted 
particular attention due to the accessibility of its Dirac 
point within a relatively large ^300 meV bulk band gap 
[4], and the availability of additional quantum well sur¬ 
face states [5-13] which may be spin-polarized by a large 
Rashba effect [14, 15]. However, spin-polarized quantum 
well states (QWS) depend sensitively on adsorbents, and 
have not been observed on clean surfaces. More generally, 
little is known about the empty state band structure of 
Bi 2 Se 3 , which is inaccessible to angle-resolved photoemis¬ 
sion spectroscopy (ARPES). For both topological surface 
states (TSS) and QWS, it remains crucial to character¬ 
ize the high energy extent to which they remain linearly 
dispersing, spin-polarized, and protected against surface- 
bulk scattering [16, 17] and backscattering [18]. 

Scanning tunneling microscopy (STM) can provide real 
space images of both filled and empty states and their 
local relationship to surface and near-surface impuri¬ 
ties. STM also provides access to momentum space in¬ 
formation via quasiparticle interference (QPI) imaging. 
When quasiparticle states of energy 6 scatter elastically 
from impurities, the interference between initial and fi¬ 
nal quasiparticle wavevectors ki and kf can result in a 
standing wave pattern with wavevector q = kf — ki at 
energy 6. The observed dispersion of q{e ), and the in¬ 
version of q{e) to find k(s), has established QPI as a re¬ 
liable k -space probe [19]. Indeed, QPI imaging has been 
used to demonstrate the protection against backscatter¬ 


ing [20, 21] and the high energy breakdown of linear dis¬ 
persion in pristine Bi 2 Te 3 [22], as well as the onset of 
backscattering in Fe cc Bi 2 _ iC Te 3 [23]. However, dispersing 
QPI in Bi 2 Se 3 has been elusive [17, 24, 25], appearing 
only in a limited energy range with inconsistent velocity 
[26-28]. Therefore, the high energy band structure and 
scattering mechanisms of Bi 2 Se 3 are unknown. 

Here we address these issues using STM measurements 
of bulk Fe cc Bi 2 _ ;E Se 3 , coupled with density functional the¬ 
ory (DFT). Previous experiments focused on the mag¬ 
netic properties of bulk Fe dopants via transport [29], 
ARPES [30], or /iSR and magnetization measurements 
[31]; or adsorbed Fe [6, 32-34]; or Fe dopants in thin 
films [35]. In this first STM study of bulk Fe^Bi 2 _ x Se 3 , 
we report multiple dispersing QPI modes, which give the 
first access to the full high energy band structure, up 
to 1 eV above the Dirac point. Furthermore, we provide 
the first real space evidence for the coexistence of QWS 
and TSS on the pristine surface of a bulk material. Fi¬ 
nally, we demonstrate the absence of backscattering of 
the QWS, and we report the relationship of scattering 
between TSS, QWS, and bulk. All of these observations 
are enabled by interstitial Fe in the van der Waals (vdW) 
gap between adjacent Bi 2 Se 3 quintuple layers (QLs). 

We studied Fe a ,Bi 2 _ a ,Se 3 single crystals with nominal 
x = 0 [36] and 0.07 [31]. Samples were cleaved in vac¬ 
uum at T ~ 40 K and immediately inserted into our 
home-built STM at 4K. Topographic images were ob¬ 
tained in constant current mode (J s ) at a fixed sample 
bias (V s ). Differential conductance dl/dV , proportional 
to the local density of states, was measured at a fixed tip- 
sample distance using a standard lockin technique. The¬ 
oretical calculations were performed using the linearized 
augmented-plane-wave method in the WIEN2K packages 
[37] and an ab initio DFT slab method [38]. 


2 



FIG. 1. (color online), (a-c) STM Topographies of pristine Bi 2 Se 3 (a, left) and three regions of Feo.o 7 Bii. 93 Se 3 [(a, right), (b), 
and (c)]. Six types of impurities, circled in (a)-(c), are segregated into three distinct micron-scale regions in Fe-doped samples. 
Region I is characterized by the same impurities seen in pristine samples, i.e. Se vacancies (Vse) and Bis e antisites. Region II 
contains Fe substitutions at two different Bi sites (Fesii and Feei 2 ). Region III contains two types of Fe interstitials (type-A 
and type-B). (d), (e) High resolution topographies of type-A and type-B Fe interstitials. Their spatial signature has a similar 
vertex-to-vertex distance of 6ao, indicated by red dots, where ao is the Se-Se lattice constant, (f) Crystal structure of Bi 2 Se 3 
with Fe dopants shown in red. Experimental conditions: (a, left) 14 = 0.4 V, / s = 30 pA; (a, right) 14 = 0.5 V, / s = 50 pA; (b) 
V s = 1V, I s = 20pA; (c) V s = 0.55 V, I s = lOOpA; (d, e) 14= 0.5V, I s = 100pA. 


Topographic images of flat terraces of Fe^B^-^Ses re¬ 
veal multiple species of three-fold symmetric defects, seg¬ 
regated into three distinct micron-scale regions, and typ¬ 
ified in Figs, l(a-c). Region I is dominated by the same 
Se vacancies and Big e antisite defects seen in “pristine” 
Bi 2 Se 3 [39, 40]. Region II is dominated by Fe substi¬ 
tutions in the top two Bi layers, Fepii and Fepi 2 [35]. 
Region III is dominated by two larger defects, detailed in 
Figs, l(d-e). The first (“type-A”), centered between top¬ 
most Se atoms, was identified as an interstitial Fe [35]. 
The second (“type-B”), centered on a topmost Se atom, 
has not been previously observed. Based on their vertex- 
to-vertex distance of 6ao = 2.4 nm, and the assumption 
that their perturbations propagate primarily along the 
ppcr chains extending out and upwards to the surface [39] , 
we conclude that both defects in region III are Fe inter¬ 
stitials at inequivalent sites within the vdW gap beneath 
the top QL. The density of atomic Fe defects in regions 
II and III is lower than the nominal doping by a factor 
of ^ 100, but occasional Fe clusters do appear [41]. 

The low Fe concentration is consistent with the small 
size of Fe 3+ compared to Bi 3+ , which makes it difficult 
to dope [42]. The spatial separation of distinct impu¬ 
rities can be understood from calculations showing that 
Fe substitution in Bi 2 Se 3 is allowed only in Se-rich con¬ 
ditions [42, 43], which also inhibit the formation of Se 
vacancies [44, 45]. On the other hand, Fe interstitials are 
not energetically favored under any Bi/Se ratio. Such in¬ 
terstitials may appear due to kinetic constraint where the 
local growth temperature is too low for Fe to overcome 
the energy barrier to substitute for Bi [42]. 

We now focus on the electronic structure in region III. 
The dl/dV maps in Figs. 2 (a-c) display clear QPI pat¬ 
terns, particularly around the new type-B interstitials. 
The QPI has three remarkable features. First, the scat¬ 
tering modes appear at unexpected high energy levels, 


almost 1 eV above the Dirac point. At such high ener¬ 
gies the band structure of Bi 2 Se 3 was expected to con¬ 
sist of bulk continuum states and was not expected to 
have a well-defined TSS [20, 46]. Second, the Fourier 
transform dl/dV maps in Figs. 2(d-f) show no dispersing 
Y-K modes, but two Y-M modes (qi and g 2 ) which dis¬ 
perse with similar slope but ^0.1 eV relative offset. To 
better isolate the dispersing Y-M modes from the long- 
wavelength dopant disorder (presumed to be isotropic), 
we plot I(qr-M ) — I{qr-K ) as a function of bias voltage 
V s and wavevector q in Fig. 2(g). Third, the Y-M mode 
velocities, v q M ~ 1.1 —1.3 eVA, are substantially smaller 
than the ARPES-measured TSS velocity near the Dirac 
point (udp = 3.5 eV A) [47]. Finally, we emphasize that 
QPI is observed only in region III. 

We argue that the parallel Y-M scattering modes in 
Fig. 2(g) demonstrate a remarkable extension of both 
TSS and QWS, far above the energy range accessible to 
ARPES [5-15]. To support this claim, Fig. 3(a) shows our 
calculated band structure in a 5QL slab, the minimal sys¬ 
tem which prevents interactions between top and bottom 
surfaces, but also supports QWS with depth comparable 
to fits of ARPES data [5, 7, 14] and the independently 
measured bulk screening length [48]. (We also checked 
that the calculated band velocity here is a robust feature, 
nearly independent of the slab thickness [41].) Constant 
energy contours (CECs) in Figs. 3(c-e) show the hexag¬ 
onal warping [49] which emerges at high energy in both 
TSS and QWS, along with the expected spin texture of 
these states [15]. Similar to previous studies on Bi 2 Te 3 
[20-22], scattering in Bi 2 Se 3 is expected to depend on the 
warping of the CECs, with a dominant mode emerging 
between two adjacent M-oriented corners of the hexa¬ 
gram (denoted by ^m)- Because these corners are sepa¬ 
rated by an angle of 7r/3, the velocity of qM should well 
approximate the band velocity along the Y-M direction: 
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FIG. 2. (color online), (a-c) dl/dV maps over a 30 x 30 nm 2 region at the indicated bias, (d-f) Fourier transform (FT) dl/dV 
maps at varying bias over a 60 x 60 nm 2 region encompassing the area of (a)-(c). A white arrow points to the dispersive T-M 
scattering mode. Inset in (d) is a schematic of the surface Brillouin zone with the high symmetry points. Data was acquired 
at T = 6 K, V s = 0.6 V, I s — 300 pA, F rms = 7 mV (a-c), V rms = 5 mV (d-f). (g) Fourier linecut showing q\ and q<i scattering 
modes. The Fourier intensity of the linecut in the T-K direction is subtracted from the Fourier intensity of the linecut in the 
T-M direction for the FT -dl/dV maps shown in (d)-(f). The white-dashed lines highlight the linear dispersions of qi and q\ 2 . 


vq M — f/c(r-M) [41]. Furthermore, the warping of TSS 
and QWS in Figs. 3(c-e) evolves similarly, but shifted in 
energy, akin to the parallel scattering modes observed in 
Fig. 2(g). Indeed, Fig. 3(b) shows excellent quantitative 
agreement between our STM scattering data and DFT 
calculations, confirming that q\ and q \2 modes arise from 
scattering between adjacent M-oriented corners of the 
CECs of the TSS and QWS1, respectively. 

Spectroscopic evidence further supports the presence 
of QWS. In Fig. 4(a), the spatially averaged dl/dV spec¬ 
trum from region III shows three distinct kinks over the 
same energy range as the observed QPI. These kinks are 
qualitatively well reproduced in Fig. 4(b), which shows 
the momentum-integrated density of states from the 
DFT band structure in Fig. 3(a). By comparison with 
Fig. 3(a), we see that the three kinks in the DOS occur 
near the three energies where the linear dispersions of the 
TSS, QWS1, and QWS2 bands break down. 

Despite the breakdown of linear dispersion and the on¬ 
set of bulk bands, both TSS and QWS remain unexpect¬ 
edly robust against backscattering and surface-bulk scat¬ 
tering. We observe no dispersing T-K modes throughout 
the measured energy range, which rules out backscatter¬ 
ing and indicates the spin polarization of the bands to al¬ 
most 1 eV above the Dirac point. We observe one weak, 
non-dispersing T-K mode around V s = 300 — 400 meV 
above [41]. Previous studies ascribed non-dispersing QPI 
modes to surface-bulk scattering in Bi 2 Se 3 [17]. However, 
our observed wavevector qx ^ 0.44 A corresponds very 
nearly to the 2\/3 a distance from the bright center to 
the lobes of the type-B interstitial Fe defects, which is 


determined simply by the bond geometry of the impu¬ 
rity [39] . We therefore find no evidence of backscattering 
or surface-bulk scattering deteriorating the TSS or the 
QWS at high energy in Bi 2 Se 3 . 

Finally, we comment on the origin of QWS. Previously 
observed QWS in Bi 2 Se 3 were believed to arise from band 
bending due to charged adsorbents [7-9, 14, 15], or vdW 
gap expansion due to burrowing adsorbents [10, 50]. Here 
we have reported the first evidence of QWS at the clean 
surface of a bulk crystal. Fe interstitials may induce 
QWS at the clean surface of bulk Fe a ,Bi 2 _ cc Se 3 by two 
mechanisms: (1) their ionization to Fe 3+ would donate 
three electrons to induce downward band bending [35]; 
(2) their cores would physically expand the vdW gap. Ei¬ 
ther mechanism supposes vertical dopant inhomogeneity, 
in keeping with the lateral inhomogeneity demonstrated 
in Fig. 1. (The tip itself can’t be responsible for the the 
observed QWS, because the tip’s positive bias can’t cause 
the downward band bending required to induce QWS in 
the bulk conduction band.) 

In conclusion, our STM/STS study of Fe^B^-^Ses 
provides fundamental new information about the empty 
state band structure of Bi 2 Se 3 . Scattering from a previ¬ 
ously unobserved interstitial Fe defect allowed QPI imag¬ 
ing and band structure reconstruction up to 1 eV above 
the Dirac point, from which we draw three main conclu¬ 
sions. First, quantitative comparison between our QPI 
data and DFT calculations proves the unexpected ap¬ 
pearance of QWS, coexisting with TSS, at the clean sur¬ 
face of a bulk material. These well-defined surface states 
persist to almost 1 eV above the Dirac point, overlap- 
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FIG. 3. (color online), (a) Calculated band structure of 
a 5QL slab with the TSS and first two QWS in blue, green 
and red, respectively, (b) Comparison between the calcu¬ 
lated high-energy Y-M bands and the measured dispersions 
of the two QPI modes. The grey circles and horizontal bars 
represent the centers and widths of Gaussian fits to the in¬ 
tensity profile at each bias in Fig. 2(g). The measured Dirac 
point bias Vdp is -0.35 V [41]. The energy offset between left 
and right axes represents a rigid band shift which arises from 
our approximation of the near-surface potential gradient as a 
square well in our DFT calculations, (c-e) Constant energy 
contours at energies spanning the transition from dominant 
TSS scattering (qi) to dominant QWS1 scattering (g 2 ). The 
color scale is defined by the density of states of the top five 
atomic layers in a 5QL slab. The expected spin texture due 
to Rashba-like spin splitting (not included in the calculation) 
is sketched by white arrows. 




FIG. 4. (color online), (a) Spatially averaged dl/dV spec¬ 
trum (U S =0.7V, I s = 0.34 nA, and V rms = 5.6 mV) showing 
three prominent kinks at 0.32, 0.45 and 0.56 V. The inset 
shows the spectrum after subtracting a third order polynomial 
background, (b) Calculated density of states (DOS) from the 
bulk (dashed line) and a 5QL slab (solid line) of Bi 2 Se 3 (with 
Vdp = 0 in the calculation). The slab DOS profile has kinks 
similar to those observed in the experimental spectrum in (a). 


ping by rsj 0.6 eV with the bulk conduction band, but 
showing no evidence of backscattering. Second, distinc¬ 
tive kinks in dl/dV spectra reveal the breakdown of lin¬ 
ear dispersion for TSS and the first two QWS. The new 
high-energy velocity is 3 x smaller than the velocity near 
the Dirac point, and is consistent with our DFT calcu¬ 
lations. Third, the absence of dispersing Y-K scatter¬ 
ing demonstrates the protection of both TSS and QWS 
against backscattering, well beyond the breakdown of lin¬ 
ear dispersion, even in the presence of sparse Fe dopants. 
Our results suggest that interstitial defects in the vdW 
gap of Bi 2 Se 3 can lead to the coexistence of protected 
spin-polarized QWS and TSS over a wide energy range, 
far outside the bulk band gap. These discoveries bode 
well for the controlled use of spins in future devices. 
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FIG. SI. An example of occasionally observed Fe clustering in Fe cc Bi 2 - ;r Se 3 crystals. The 30 x 30 nm 2 topography was acquired 
at T = 7K, V s = 0.3V, I s = 30pA. 
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FIG. S2. (a) Topography of a 30 X 30 nm 2 region of Fe^Bh-xSes (T = 6K, V s = 0.6 V, I s = 30 pA). (b-f) dl/dV maps 

over the region shown in (a) at the indicated bias. Interference waves are seen around the type-B defects. A white arrow is 
overlaid on the five images and shows that as the energy is increased the wavefront moves towards the impurity (T = 6K, 
V s — 0.6 V, I s — 300 pA, V rms = 7mV), (g) Schematic diagram of the Brillouin zone (black) and a hexagonally warped constant 
energy contour appropriate for Bi 2 Se 3 . The magnitude of scattering vector qM should be similar to the value of initial and 
final quasiparticle wavevectors ki and kf because of the 7r/3 angle separation. Thus, the velocity of qm is expected to be equal 
to the Y-M band velocity, (h-1) Fourier transforms of the dl/dV maps at varying bias over 60 x 60 nm 2 regions encompassing 
the area of (a-f) (T = 6K, V s = 0.6 V, I s — 300pA, V rms = 5mV). A white arrow is overlaid to show that as energy increases 
the scattering wavevector increases in the Y-M direction. 



FIG. S3, (a, b) Linecuts of Fourier transform dl/dV spectroscopy maps shown in Fig. 2 of the main text and Fig. S2. Linecuts 

are taken in the (a) Y-M and (b) Y-K directions. The difference of the two images gives the result shown in Fig. 2(g) of the 
main text, which highlights the q\ and Y-M scattering modes. The dl/dV maps shown in Fig. S2 and Fig. 2 of the main 
text were normalized by energy layer, Fourier transformed, six-fold symmetrized, and mirrored along a Bragg direction. 
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FIG. S4. dl/dV spectrum acquired in region III with V s =0.5 V and I s = 100 pA shows a Dirac point at —0.35 V. 



Bias (V) 

FIG. S5. dl/dV spectra were acquired at the listed experimental conditions, with varying junction resistance R — V s /I s - The 
junction resistance is exponentially dependent on the tip-sample separation. Blue, green, and red vertical bars are overlaid to 
show the approximate energies of kinks in the spectra, which we identify as the breakdowns of linear dispersion of the TSS, 
QWS1, and QWS2, respectively. The kink energies are consistent across spectra taken with different experimental conditions, 
both in and out of magnetic field. Experimental conditions not listed in the legend: (a) V s — 400 mV, V rms = 14 mV; (b) 
V s — 500 mV, Vrms — 14 mV; (c) V s — 600 mV, V rms — 5 mV; (d) V s — 700 mV, V rms — 5.6 mV. Spectrum (d) is the same 
data presented in Fig. 4 of the main text. 
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FIG. S6. (a-c) DFT slab calculations with various QL thicknesses, (d-f) Band structure along the T-M direction at the 

indicated QL slab thickness, in an energy and fc-space range corresponding to the qi and q 2 scattering modes in the data. The 
black line shows the best linear fit (E — vk + Eq) in the range k — 0.2 A -1 to k — 0.3 A -1 . The velocity extracted from the 
fit is equal to ~ 1.1 eVA and is nearly independent of the slab thickness, (g-i) Calculated DOS at varying QL thicknesses. 
Dashed vertical lines in (g)-(i) are guides illustrating the energies of the kinks in the DOS. These kink energies are also marked 
as horizontal lines in (a)-(c). Although a full Poisson-Schrodinger calculation of the potential gradient profile [5, 14] is beyond 
the scope of this work, the calculated DOS kinks arising from the simple square well potential with thickness 5QL, 6QL, or 
8QL all show excellent qualitative agreement with the measured data in main text Fig. 4a. 
































